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Abstract 

We investigate the effect of supersymmetric CP violating phases on 
the B — > X s 7 and B — > X s l + l~ decays in the minimal supergravity 
model. We show that the phase of the trilinear scalar coupling constant 
for top squarks is strongly suppressed and aligned to that of the gaug- 
ino mass due to a renormalization effect from the Planck scale to the 
electroweak scale. As a result, the effect of supersymmetric CP violat- 
ing phases on the B — > X s 7 and B — > X s l + l~ decays are small taking 
into account the neutron and the electron electric dipole moment con- 
straints. For the B — > decay, the amplitude has almost no new CP 
violating phase and the direct CP asymmetry is less than 2 %. For the 
B — ► X s l + l~ decay, the branching ratio can be sizably different from that 
in the standard model only when the sign of the B — > X s 7 amplitude is 
opposite to that in the standard model. 
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The origin of the CP violation is one of main issues in current particle physics. 
In the standard model (SM) the CP violation is originated from the phase of the 
Kobayashi-Maskawa matrix [|TJ . A new source of CP violation can appear in models 
beyond the SM. 

Among various models beyond the SM, the minimal supersymmetric SM (MSSM) 
is one of the most promising candidate. The MSSM contains many new parameters, 
i.e. soft supersymmetry (SUSY) breaking parameters. If we allow arbitrary soft 
SUSY breaking parameters, too large flavor changing neutral current (FCNC) pro- 
cesses, such as K°K° mixing, appear. In the minimal supergravity (mSUGRA), soft 
SUSY breaking parameters are assumed to have universal structures at the Planck 
scale, so that these dangerous FCNC processes are suppressed. 

In the mSUGRA, there is no intrinsic reason that these SUSY breaking param- 
eters should be real and there are four new CP violating phases, i. e. phases of the 
gaugino mass, the higgsino mass parameter, the SUSY breaking Higgs boson mass, 
and the trilinear scalar coupling constant, of which two combinations are physically 
independent. These phases induce the neutron and electron electric dipole moments 
(EDMs). There are many works on the constraints of the EDMs in the MSSM as 
well as in the mSUGRA It is shown that in the mSUGRA, if we take a phase 

convention that the trilinear scalar coupling constant and the higgsino mass param- 
eter have phases, 4>a and </> M , respectively, the constraint on cj)^ is much stronger 
than that on (ft a [@HI • 

It is known that the B — > Xg'-f process gives strong constraint on the SUSY 
model. In particular, in Ref. J7|, rare B decays, such as B — > Xgj and B — ► 
X s l + l~, are studied in the mSUGRA without new CP violating phases. For the 
B — > Xg'y decay, the SUSY contributions interfere with the amplitude in the SM 
either constructively or destructively, and the amplitude can change its sign. It is 
also shown that the B — > X s l + l~ branching ratio is enhanced compared to the SM 
prediction if the sign of the B — > Xg'j amplitude is opposite to that in the SM. It is 
interesting to investigate the effect of the CP violating phases on various B decays. 
In Ref. one of the authors (T.N.) analyzed effect of the SUSY CP violating phase 
on B°B° mixing, and showed that the effect is small. Recently a possibility of large 
direct CP asymmetry in the B — > X s ^ process is studied in the MSSM [0,0], and 
MSSM with SUGRA-motivated SUSY breaking terms P . 
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In this letter we investigate the effect of the SUSY CP violating parameter 
on rare B decays, B — > X s ^ and B — > X s l + l~, in the mSUGRA. In our analysis 
we require the universality of SUSY breaking terms at GUT scale and explicitly 
solve the renormalization group equations (RGEs) to determine the masses and the 
mixings of SUSY particles and also require the condition for the radiative electroweak 
symmetry breaking. We confirm that M is strongly constrained by the neutron and 
electron EDM bounds whereas 4>a is almost unconstrained. However we show that 
the phase of the A-term for top squarks is reduced due to the large top Yukawa 
coupling constant and aligned to that of the gaugino mass. Therefore the phase of 
the A-teim. for top squarks is strongly suppressed. We show that the CP asymmetry 
in rare B decays is suppressed in the mSUGRA if the neutron and electron EDM 
constraints are taken into account. 

In the MSSM, the Yukawa coupling constants are described by the following 
superpotential, 

IUmssm = {YuhQiUjHi + (Y D ) ij Q i D j H 1 + {Ye^E.L^ - fiH^, (1) 

where Q and L denote the SU(2) L quark and lepton doublets, U, D, and E are 
SU(2) L singlets, and Hi, H 2 are SU(2) L Higgs doublets. The i,j represent genera- 
tion indices. In addition to the SUSY invariant terms, there are following soft SUSY 
breaking terms, 

- £ soft = (m^ijqliqLj + (mDijU^URj + (m 2 D ) i:j d* m d Rj 

+ ( m l)ij Ai 4? + ( m Eh e Rj 

+ Alhlh + A 2 2 hlh 2 + {Biih x h 2 + H.c.) 

+ [(A^ijqL^ u* Rj + {A^ijhiqLid*^ + (A E )ijhie Ri l L j + H.c. 

+ (-M 1 BB + -M 2 WW + -M 3 GG + R.c). (2) 
2 2 2 

Hereafter we denote the superpartners by small letters with tilde. The SUSY break- 
ing terms depend on details of SUSY breaking mechanism. In the minimal su- 
pergravity model these SUSY breaking terms are originated from the gravitational 
interaction and given by the following universal structure at the high energy scale: 

M 1 = M 2 = M 3 = M x , (3) 
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( m Q)ij = ( m u)ij = ( m D)ij = ( m L)ij = ( m Eh = m S ij: (4) 

A? = Aj = ml (5) 

(Au)ij = A x (Yu) ij: (A D )ij = A x (Y D )ij, (A E ) lj = A x {Y E ) ij . (6) 

For simplicity we assume the GUT relation for gaugino masses and put the universal 
condition at the GUT scale (~ 2 x 10 16 ) neglecting the renormalization effect from 
the Planck scale to the GUT scale. The SUSY breaking parameters at the elec- 
troweak scale are obtained by solving the RGEs. In principle, the parameters, M x , 
A X) /i, and Bfi, can have phases. Since only two combinations of the four phases 
are physical CP violating phases, we take only A x and /i as complex parameters 
hereafter. 

In order to see qualitative feature of RGEs for A-terms let us first neglect flavor 
mixings in the RGEs. The RGEs of A-terms for the first and second generations are 
given by 

^-A ei = 3a 2 M 2 + 3a 1 M 1 - a T A T - 3a b A b , (7) 
at 

-j;A di = ^-a 3 M 3 + 3a 2 M 2 + - a T A T - 3a b A b , (8) 

at 6 9 

d 16 13 

—A u . = —a 3 M 3 + 3a 2 M 2 + —a 1 M 1 -3atA t , (9) 
where i — 1,2 and for the third generation 

~j~A T = 3a 2 M 2 + 3a 1 M 1 - Aa T A T - 3a b A b , (10) 
d 16 7 

— A b = —a 3 M 3 + 3a 2 M 2 + -a>iM 1 — a T A T — 6a b A b — a t A t , (11) 

4 A = ^a 3 M 3 + 3a 2 M 2 + ^-a l M l -a h A b -Qa t A t . (12) 
at 3 9 

Here A^ = (Af)a/ (Yf)ni t = — ln(Q 2 ) / (An) where Q is a renormalization point, = 
gf/(Aix), and a ft — Y^./(47r). In the right-hand sides (RHS's) of above equations, 
only the Yukawa coupling constants of the third generation are retained. Since the 
RHS's of RGEs for A-terms depend linearly on the A-terms and gaugino masses, 
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general solutions can be written in terms of the universal A-term and gaugino mass 
as follows: 



A fi =CfA x -C 9 fi M x , (13) 

where the coefficients Cf. and Cj. are functions of the Yukawa and gauge coupling 
constants. In Fig. [I], Cj. and Cf. are shown as a function of tan/3 (= v 2 /v 1 ). Cf 
is much smaller than C 9 because Cf is reduced by the large top Yukawa coupling 
constant. Therefore, the phase of A-term for top squarks is strongly suppressed due 
to the renormalization effects even if the phase of A-term at Mx scale is maximal. 
Considering the current experimental lower bound on the chargino mass, m%+ > 91 



GeV ||11|| , Mx at the GUT scale must be roughly larger than 120 GeV. In principle, 
the contribution from Ax can dominate in Eq. ( |T3"D if Ax is larger than O(10) 
TeV. However, it makes scalar particles heavier than 1 TeV, in which case SUSY 
loop effects on FCNC processes in B decays are small. As tan/3 becomes larger, 
Cf. for the bottom squark is reduced due to the bottom Yukawa coupling constant. 
Note that the suppression of the phase of A t is general feature in models where the 
A-terms are generated at a high energy scale. 

Let us discuss phenomenological consequences of SUSY CP phases on the 
B — > Xg'-f and B — > X s l + l~ decays. These processes are described by the following 
effective Lagrangian 

AC 10 

C eff = —£v tb V t * s Ci(Q)Oi(Q). (14) 
V2 i=1 

The operators, O 7 -Oi , are most relevant for the calculation of the processes, which 
are given by, 

°7 = -r^rn b sIa^b R F^ (15) 
Os = -^m b !Za^T a b R F« u , (16) 

9 = rf^slTAW (17) 

l07T z 

Oio = ^-^llM^ld- (18) 
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In order to calculate the B —>■ X s j and B —>■ X s l + l~ processes in the mSUGRA, we 
first solve the RGEs of the MSSM to determine the masses and mixings of SUSY 
particles. Then, integrating out SUSY particles at the electroweak scale, the SUSY 
contributions are included into the Wilson coefficients Cj in matching conditions. 
The Wilson coefficients at the bottom mass scale are calculated by solving the RGE 
of QCD at the next-to-leading order (NLO). As for the NLO calculation we follow 
the results in Ref. 



121 for the B 



B 



X s 7 process, and the results in Ref. |T3[ for 
XJ + l~ process. 

The direct CP asymmetry in the B 



A CP (S) 



T(B -> X sl ) - T{B 



X s 7 decay is given by 



T(B - 

Q3 (/■*&) 
|CV| 2 



X al ) + T(B^Xrf) 



40 



-[v(z)+b(z,5)]Im[(l + e s )C 2 C* 7 ] 
y 



--Im(C 8 CV*) + -b(z, S)lm [(1 + e s )C 2 C* 8 ] 



(19) 



where 5 is an energy cutoff parameter for the photon, fi b is a renormalization point 
at the bottom mass scale, z = (m c /m b ) 2 , e s = V ub V* s / (V tb V t * s ) , and functions v and 
b are found in Ref. ||. In the SM, the CP asymmetry is estimated as A^p{5 = 
0.99) ~ 1.5 X 10 2 ?7 where r\ is the Wolfenstein parameter. The SM prediction is 
small because the small parameter e s , which is O(10 -2 ), is the only source of the 
direct CP violation in the B — ► X s ^ process. If C 7 or C 8 has a sizable new CP 
violating phase, Aqp could be large. 

The dilepton spectrum of the B — > X s l + l~ decay can be written by 



dB(B 



XJ+l 



ds 



B(B -> XJv) 



2 v tb v; s * 



Att 2 



1 



X 



e//|2 i \ r i I 
+ WO 



! )(l + 2j) 



V cb f(m c /m b )K(m c /m b ) 
2 (2 + s) + 12Re(CV*C 



(20) 



where s is the dilepton invariant mass square normalized by bottom mass square, 
/ = l - Sx 2 + 8x 6 - x 8 



24x In x, and k is a QCD correction factor [pL3| . Since 



there is an interference term between CV and , the dilepton spectrum depends 
on the phases of Cj or Cg . 
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For numerical calculations in the mSUGRA, we scanned the SUSY parameters 
in the range < tuq < 1 TeV, < Mx < 0.5 TeV, \Ax\ < 5mo, and we follow 
Ref. |7[| for detailed procedures of the calculation. For definiteness, we take the 
KM parameters, |V US | = 0.2196, \V cb \ = 0.0395, \V ub /V cb \ = 0.08, and <J 13 = vr/2 in 
the standard parametrization [14|j. As for the B — > X s 7 decay, we take 5 = 0.99, 
fib = m b = 4.8 GeV, and m t = 175 GeV. 

In Fig. H the neutron EDM is shown as a function of 0^ and 4>a for tan f3 — 30. 
In the numerical calculation only the neutron EDM from the quark EDMs is in- 
cluded. Recently it was pointed out that the EDM constraints may be relaxed by a 
cancellation among different contributions in the mSUGRA However, we do 

not rely on such a cancellation because each contribution has different hadronic un- 
certainty, so that it is difficult to determine the parameters where such a cancellation 
occurs. We are only interested in the case where SUSY particles are lighter than 1 
TeV because otherwise SUSY effects on the B — > and B — > X s l + l~ decays are 
strongly suppressed. In such pointed out in Ref. f§-§|, the neutron EDM 

exceeds the present experimental bound, \d n \ < 0.97 x 1CT 25 e cm |15[], unless 0^ is 



< 10 2 . On the other hand, 4>a can be 0(1) even if the masses of SUSY particles 
are O(100) GeV. 

The real and imaginary parts of CV at the bottom mass scale divided by the 
SM value are plotted in Fig. § for tan/3=3, 10, 30. As in the case of no SUSY CP 
violating phase, the SUSY contributions to C 7 and C s become large, however, those 
to Cg and C*io are small. In this figure, the experimental bound on the B — > X s ^ 
branching ratio, 2.0 x 10~ 4 < B(5 — > X^) < 4.5 x 10~ 4 |16 l, is imposed, therefore 
the region between two circles are allowed. It is interesting to see that without the 
neutron and electron EDM constraints, C7 can have different phase from the SM 
value. On the other hand, with the EDM constraints, i.e., \d n \ < 0.97 x 10~ 25 e 



cm, \d e \ < 4.0 x 10 e cm IpTj , the imaginary part of Cj/Cj is quite suppressed 



and either CV ~ Cf M or Cj ~ — Cf M region is allowed. This is a similar result 
to that obtained without the SUSY CP violating phases [0. It is known that the 
charged Higgs boson and the chargino contributions to CV can be significant, and 
that the charged Higgs boson contribution to C7 has the same phase as the SM 
contribution. On the other hand the chargino-stop loop contribution to C7 depends 
on the new SUSY CP phases. In order to have the large phase of CV, the imaginary 
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part of the chargino contribution must be large. With the neutron and electron 
EDM constraints, however, M must be quite small. Moreover the phase of A t is 
also suppressed due to the RGE effect as discussed above. Therefore it is difficult to 
have large phase of C-j. We also find that Cg/Cg M does not induce large imaginary 
part although the magnitude itself can be changed by the SUSY contributions. 
This means that large CP violating phases do not appear in rare B decays in the 
framework of the mSUGRA even if the new CP violating phases are introduced. 



This is a distinct feature from the result which is obtained in Ref. |Ij where the 
authors did not follow the universality condition for scalar masses at the GUT scale. 
From Fig. 3(a)-(c), only for tan/3 =30, there is parameter space where C-jjC^ M is 
negative. We find that in this parameter region the lighter stop mass is less than 
about 200 GeV as in the case of no new SUSY CP phase J7|. 

In Fig. ^(a)-(b), we plot the B — > X s fi + fi~ branching ratio and the direct CP 
asymmetry in the B — > X s 7 versus the neutron EDM for tan/5 = 30. After taking 
into account the neutron and electron EDM constraints, we show that there are 
two branches of the B — > X s [i + fi~ branching ratio and the larger branching ratio 
corresponds to the case where the sign of Cj is opposite to that in the SM. The 
branching ratio can be about twice as large as that in the SM with Cj ~ —C^ M . 
The CP asymmetry turns out to be less than 2 %. 

In conclusion, we investigate the effect of the SUSY CP violating phases 
(0 M , on the rare B decays in the mSUGRA model taking into account the RGEs 
for the SUSY breaking parameters. If the SUSY particles are in the hundred GeV 
region, M is strongly constrained by the EDM bounds. On the other hand, the 
phase of A-term for top squarks is aligned to that of the gaugino masses due to the 
RGEs. As a consequence, the effect of the SUSY CP violating phases is small and 
either C-j ~ Cf M or C-j ~ — Cf M is allowed. We show that the direct CP asymmetry 
is less than 2 % taking into account the EDM constraints. For the B — > X s l + l~ 
decay, there is a twofold ambiguity of the branching ratio according to the sign of 
C-j. The branching ratio can be twice as large as the SM value when Cj ~ — Cj M . 
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Figure Captions 



FIG. 1 The coefficients C A and C 9 in Eq. (0) are shown as a function of tan f3. 
Here we take MS masses for quarks as m u = 3.3 MeV, = 6.0 MeV, and 
m s = 120 MeV at the scale of 2 GeV. We take pole masses as m c = 1.4 GeV 
and m b = 4.8 GeV, m t = 175 GeV. 

FIG. 2 The absolute value of the neutron EDM (\d n \) is plotted as a function of 
M (a) and <$a (b) for tan f3 = 30. Here input SUSY parameters are scanned in 
a region, < m < 1 TeV, < Mx < 0.5 TeV, and \Ax\ < 5m . The dashed 
line represents the present experimental upper bound, \d n \ < 0.97 x 10~ 25 e 
cm. For Fig. 2(b), squares correspond to the parameter spaces M = 0,tt. 

FIG. 3 Cf/Cj M at the bottom mass scale is shown imposing the current exper- 
imental bound for the B —>■ Xgj branching ratio for tan/3 = 3 (a), 10 (b), 
30 (c). Dots correspond to values without the neutron and electron EDM 
constraints and squares correspond to values with the EDM constraints. The 
input SUSY parameters are the same as Fig. 2. 

FIG. 4 (a). The B — > X s fi + fi~ branching ratio is plotted as a function of the neu- 
tron EDM. As for the B — > X s fi + fi~~ branching ratio, in order to avoid the 
J/ip resonance, we integrate the dilepton spectrum in a region 4m 2 < s < 
( m J/ip — O.l(GeV)) 2 where s is the dilepton invariant mass square, (b). The 
absolute value of the direct CP asymmetry in the B — > X s , y (\Acp\) is plotted 
as a function of the neutron EDM. In these figures the input SUSY parameters 
are the same as Fig. |2|. The dashed line represents the present upper bound 
for the neutron EDM. Squares correspond to the parameter spaces with the 
electron and neutron EDM constraints. 
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